Electronic structure and magnetic anisotropy of Cr02 
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The problem of importance of strong correlations for the electronic structure, transport and mag- 
netic properties of half-metallic ferromagnetic Cr02 is addressed by performing density functional 
electronic structure calculations in the local spin density approximation (LSDA) as well as using 
the LSDA+U method. It is shown that the corresponding low-temperature experimental data are 
best fitted without accounting for the Ifubbard U corrections. We conclude that the ordered phase 
of Cr02 is weakly correlated. 



PACS numbers: 71.27.+a 75.30.Gw 79.60.-i 

As a compound with multiple industrial applications 
and its unusual half-metallic electronic structure, Cr02 
has recently attracted a lot of theoretica l^i^i^i^i^i^ and 
experimental2i2i2iiSiii interest. The main discussion was 
centered around the role of strong correlations for the de- 
scription of its ferromagnetic phase. Since Cr in its for- 
mal 4-f valence state has two id electrons of t2g symme- 
try, one would expect manifestation of correlation effects 
of the Mott-Hubbard nature. On the other hand, metal- 
lic behavior of spin majority band suggests that Coulomb 
interactions of the Hubbard type can be screened out^. 
The comparison with the available photoemission and 
optical conductivity data did not make the situation 
more clear. One-electron spectra calculated using the 
LSDA-I-U methodiSii^ fit well the photoemission and in- 
verse photoemission experiments with the choice of intra- 
atomic Coulomb and exchange parameters ?7 = 3 eV and 
J = 0.87 eV^. This indicates the importance of strong 
correlations. Contrary to this result, the LSDA opti- 
cal conductivity calculations explain experimental data"*, 
which suggests the regime of weak coupling. 

In the present paper we address the issue of contro- 
versial role of strong correlations in ferromagnetic Cr02 
by presenting combined studies of its electronic struc- 
ture, optical conductivity and magnetic anisotropy us- 
ing the LSDA and LSDA-I-U schemes. We employ a 
lincar-muffin-tin-orbital (LMTO) method in its atomic 
sphere approximation (ASA)^^*^^ for our electronic struc- 
ture calculations. The low symmetry of the rutile struc- 
ture and small packing factor of the unit cell require an 
introduction of additional empty spheres. Their posi- 
tions are chosen to be 4c and Ag in Wyckoff notations. 
The radii of the spheres (in atomic units) for Cr and O 
atoms, as well as of the empty spheres are chosen to be 
1.975, 1.615, 1.378 and 1.434, correspondingly The ba- 
sis set adopted in the calculations is Cr(4s, 4p, 3d) and 
0(2s, 2p). 



In rutile structure Cr atoms are surrounded by dis- 
torted oxygen octahedra. The positions of the octahedra 
lead to a new natural basis for Cr orbitals. In this ba- 
sis the cubic component of the octahedral crystal field 
splits the fivefold degenerate 3(i orbital into higher en- 
ergy doubly degenerate eg level and lower energy triple 
degenerate t2g level. Distortions of oxygen octahedra fur- 
ther split the t2g states into lower energy t\g orbital [xy 
character) and higher energy twofold degenerate or- 
bitals {yz -\- zx and yz — zx characters)^. 

The results of the LSDA band structure calculation in 
the vicinity of the Fermi energy are shown in Figs. ^ and 
13 The Fermi level crosses the spin majority t2g mani- 
fold. The rest of the Cr 3d-states is formed from four 
eg-bands and three t2g spin minority bands which are 
located above the Fermi level. In both spin channels Cg 
and t2g bands are well separated for all momenta ex- 
cept for the F-point. The whole 3fi-complex is strongly 
hybridized with oxygen. In Fig. [5] one can see that in 
the spin minority channel there is gap of approximately 
1.3 eV between the oxygen 2p-band and the chromium 
d-band. This gap leads to 100% spin polarization at Ep 
and assures the magnetic moment to be precisely equal 
to 4 /is per unit cell. The t2g bands that cross the Fermi 
level in the spin majority channel mainly consist of the 
orbitals (see Fig. Almost non-dispersive narrow 
band below Ep (shown as lightly shaded) is formed by 

the t\g orbital. This localized state undergoes large ex- 
change splitting A*^^ making spin minority i|g orbitals 
unoccupied (see Fig. I^J. 

The main changes which occur in the band struc- 
ture for non-zero values of U and J using the LSDA-I-U 
method are schematically shown in Fig. O These calcu- 
lations were performed with [/ = 3 eV and J = 0.87 eV. 
The center of gravity of occupied t\g band is pushed down 
by 0.6 eV. The spin minority unoccupied Cg bands are 
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FIG. 1: LSDA band structure of Cr02 for spin majority 
carriers. Dark and light shaded areas show the specific weight 
of t^g and t|g orbitals respectively in the particular band. 
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FIG. 2: LSDA band structure of Cr02 for spin minority 
carriers. Dark and light shaded areas show the specific weight 
of and t|g orbitals respectively in the particular band. 

pushed up by 0.6 eV, which opens 0.4 eV gap between 
and Eg bands above the Fermi level. In the spin mi- 
nority channel the occupied oxygen bands are shifted up 
by 0.3 eV. The upper unoccupied t2g and Cg bands are 
shifted up by 1.1 eV. As a result, the insulating gap is 
increased and reaches the value of 2.1 eV. 

Now we compare our calculated electronic structure 
using the LSDA and the LSDA+U method with the 
available experimental data. Fig. ^ shows comparison 
of ultraviolet photoemission spectroscopy (UPS) exper- 
iments^ (photon energy hiy = 40.8 eV) with the theo- 
retical spectra which are calculated densities of states 
smeared by both Gaussian and Lorentzian broadening 
functions. The Gaussian broadening takes into account 
experimental resolution while Lorentzian takes into ac- 
count finite lifetime effects. The Gaussian broadening 
parameter is taken to be 0.4 eV. The full width at half 



2 



LDA LDA+U 




FIG. 3: Schematic density of states (DOS) of Cr02 deduced 
from the LSDA and LSDA-I-U calculations. Shaded semicir- 
cles from right and left represent the bands for spin majority 
and spin minority carriers. 

maximum (FWHM) of the Lorentzian was taken to be 
energy dependent and equal to Q.2\E — Ef\ eV. We can 
distinguish two main features in the UPS spectra: (i) 
a small hump in around —1.5 eV which arises from the 
t2g band of Cr, and (ii) a big hump around —6.0 eV 
which comes from the broad 2p oxygen band. Both fea- 
tures are fairly well described by both the LSDA and the 
LSDA-I-U calculation. The small discrepancy between 
the LSDA calculation and experiment could be referred 
to the fact that at small photon energies photoemission is 
a more surface sensitive technique. Indeed, recent PES 
studies of Vanadium oxides^^ have been found to yield 
spectra not characteristic of the bulk, but rather of sur- 
face atoms whose lower coordination number can render 
more strongly correlated surface layer. 
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FIG. 4: Comparison between theoretical densities of states 
and experimental— UPS spectra for Cr02. The theoretical 
DOS were smeared out by Gaussian and Lorentzian broad- 
ening functions to account for experimental resolution and 
lifetime effects. The secondary electron background has been 
taken into account. 

For the unoccupied states we have chosen to compare 
our results with the available x-ray absorption spectra 
(XAS)^^ rather than with the inverse photoemission as 
it had been done before-. The main reason for this is 
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that XAS is a bulk (not surface) sensitive method. The 
2p Cr XAS spectrum^ ^ is compared to our theoretical 
calculations in Fig. El To deduce theoretical spectra we 
performed both Gaussian and Lorentzian broadening of 
3(i and 4s partial DOSes. Two first peaks around 0.5 eV 
and around 1.5 eV come from the unoccupied Zd orbitals 
of chromium. The main contribution to the second peak 
comes from the t2g orbitals in the spin minority channel. 
Thus, the LSDA+U overestimates the spin minority gap 
twice as much. 
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FIG. 5: Comparison between theory and experimenli-- for Cr 
2p x-ray absorption (XAS) spectrum. To deduce theoretical 
curve from the partial Cr 3d DOS we used 0.1 eV for Gaus- 
sian FWHM. The Lorentzian FWHM was taken to be energy 
dependent and equal to 0.2|i5 — Ef\- The binding energy 
of core 2p3/2 Cr state 577 eV has been subtracted from the 
experimental spectrum. 

Below we discuss the optical conductivity of Cr02. In 
Fig. ini diagonal x-components of the optical conductiv- 
ity calculated using the LSDA and LSDA-I-U methods 
are compared with the experimental results reported by 
Basov and coworkers— {x coordinate refers to the basis of 
unit cell). The main two features of the calculated optical 
conductivity are a shoulder around 2 — 3 eV and a broad 
hump located at energies 0.2 — 1.5 eV. In both LSDA and 
LSDA+U schemes the shoulder can be identified with two 
types of transitions. First contribution arises from the 
minority spin gap transitions and the second one comes 
from transitions between the occupied t\g and unoccu- 
pied Eg bands. The hump is formed by interband tran- 
sitions within the t2g-Tiia3i\io\d and the oxygen 2p bands 
near the Fermi level in the spin majority channel. Ap- 
parently, the LSDA prediction is much closer to the ex- 
perimental curve than the LSDA-I-U one. The LSDA-f U 
calculations overestimate the minority gap, and due to 
that, the spin minority transitions occur at higher ener- 
gies. 

Results of the calculated magnetic anisotropy of Cr02 
are presented below. Magnetic anisotropy is the depen- 
dence of internal energy on the direction of spontaneous 
magnetization. The magnetic anisotropy is a relativistic 
phenomenon arising due to spin-orbit coupling, where 
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FIG. 6: Comparison of the optical conductivity of Cr02 ob- 
tained using the LSDA and LSDA-I-U methods against the 
experimental datai^. 



the spin degrees of freedom interact with the spatial 
anisotropy through the coupling to the orbital degrees of 
freedom. This induces a preferred direction of spins. Be- 
cause Cr02 is a metastable compound, which irreversibly 
decomposes at about 200 °C, all measurements for this 
material are performed on polycrystals, microgranulus 
or thin films. The first reliable result on magnetocrys- 
talline anisotropy measurements was reported in Ref. llTl 
The discovery of the atmospheric pressure chemical va- 
por decomposition (CVD) technique has allowed to grow 
high-quality films of Cr02. As a result, in recent years a 
lot of studies of magnetic properties were performed on 
epitaxial Cr02 layers deposited on single crystal (100) 
Ti02 substrate3i2ii2*2a. For thicker films (700 A-1.2/im) 
the in-plane magnetic anisotropy was observed with [001] 
and [010] easy and hard axis directions respectively. The 
magnetocrystalline anisotropy constant Ki has been re- 
ported by different groups to be 4.4 x lO^erg/cm^ — , 
2.7 X lO^erg/cm^ M. and 1.9 x lO'^erg/cm^ — . How- 
ever, these values can significantly differ from the bulk 
quantities because of a large lattice mismatch between 
Cr02 films and Ti02 substrates (till 4% ). The relax- 
ation as a function of thickness is very gradual and even 
for 1.2/im films magnetic anisotropy shows significant de- 
viation from the bulk value. 

We calculate the magnetic anisotropy energy (MAE) 
by taking the difference of the two total energies with dif- 
ferent directions of magnetization ([001], [010], [111] and 
[102]). For the momentum space integration, we follow 
the analysis given by Trygg and co-workersSi and use 
special point method^^ with a Gaussian broadeningS^ of 
15 mRy. The validity and convergence of this proce- 
dure has been tested in their work—. We used about 
1000 k-points in the irreducible Brillouin zone, while the 
convergence of MAE is tested up to 8000 k-points. 

The direction [001] was found to be easy magnetiza- 
tion axis within our LSDA calculation which is consis- 
tent with latest thin film experimentsiSiiSiSS. Numerical 
values of MAE in this case exceed the maximum experi- 
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mental value by approximately two timesSa. 

To figure out the influence of intra-atomic repulsion U 
on the magnetic anisotropy, we have performed LSDA+U 
calculations for different values of U increasing it from 
to 6 eV (J = 0.87 eV has been kept constant except for 
the LSDA U = case). The results of these calculations 
are presented in Fig. [7| MAE decreases rapidly starting 
from the LSDA value (which is approximately equal to 
68 /X eV per cell) and changes its sign around U « 0.9 eV. 
This leads to switching correct easy magnetization axis 
[001] to the wrong one, namely [102]. The biggest ex- 
perimental value of the MAE reported in the literature 
is 15.6 fj, eV per cell^'^. The calculated MAE approaches 
this value around U = 0.6 eV. This signals that corre- 
lation effects in the d-shell may be important for this 
compound although they are strongly screened out. 
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FIG. 7: The magneto-crystalline anisotropy energies for Cr02 
as functions of U. The experimental value of MAE £[010] — 
-E[001] — Ib.Gfj, eV per cell is shown by arrow. 

To conclude, we have reported the LSDA and 



LSDA-I-U calculations of electronic structure, optical 
conductivity and magnetic anisotropy of Cr02. Our 
comparisons with the experimental data revealed that 
the best match is provided by the local spin density ap- 
proximation. We explained the discrepancy between the 
LSDA and photoemission studies, discussed earlier by 
other authors'^'^, by the fact that due to small photon 
energies used in PES, it is more surface rather than bulk 
sensitive technique. We resolved this problem by show- 
ing that XAS spectrum is unambiguously described by 
the LSDA calculation. Is has been also shown that even 
intermediate values of U{ of the order of 1-2 eV) lead to 
the failure of the LSDA-I-U method to describe the mag- 
netic anisotropy and the optical conductivity of Cr02. 
Since the LSDA-I-U is not adequate for the description 
of electronic structure of Cr02 as well as of its optical 
and magnetic properties, we conclude that the ordered 
phase of Cr02 could be described as weakly correlated 
material with small values of on-site Coulomb repulsion. 
It is important to notice that, while we have found that 
the a simple one-electron picture describes well the ferro- 
magnetic phase of this material, there is a narrow band 
formed by the non-dispersive orbitals {xy character) 
which in the paramagnetic phase will be single occupied, 
due to the on-site Coulomb interactions, an effect which 
cannot be described in LDA and will require a Dynam- 
ical Mean-Field treatment for this materials as done in 
Ref. 5. The physical basis for the applicability of static 
mean-field picture in the ferromagnetic phase of this ma- 
terial, is due to the large exchange splitting which is able 
to effectively enforce the single occupancy of the or- 
bitals. 
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